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We present a low-temperature experimental test of the fluctuation theorem for electron transport 
through a double quantum dot. The rare entropy- consuming system trajectories are detected in 
the form of single charges flowing against the source-drain bias by using time-resolved charge detec- 
tion with a quantum point contact. We find that these trajectories appear with a frequency that 
agrees with the theoretical predictions even under strong nonequilibrium conditions, when the finite 
bandwidth of the charge detection is taken into account. 
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The second law of thermodynamics states that a 
macroscopic system out of thermal equilibrium will ir- 
reversibly move toward equilibrium driven by a steady 
increase of its entropy. This macroscopic irreversibility 
occurs despite the time-reversal symmetry of the under- 
lying equations of motion. Also a microscopic system will 
undergo an irreversible evolution on a long time scale, but 
over a sufhciently short observation time r, both entropy- 
producing trajectories as well as their their time-reversed 
entropy-consuming counterparts occur. It is only because 
of the statistics of these occurrences that a long-term ir- 
reversible evolution is established. This phenomenon is 
described by the fluctuation theorem [l|, |2| . 

Irrespective of the description of the trajectories be- 
ing system-specific, the fluctuation theorem (FT) relates 
the probabilities Pr{AS) for processes that change the 
entropy of the system by an amount AS* during an arbi- 
trary time T by the equation 
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where /cg is the Boltzmann constant. In a seminal work, 
Wang et al. [SJ] tested the FT by measuring force trajecto- 
ries of a micron-scale latex bead in a liquid. Since then, 
the FT has been tested in other systems [j-lSJ, but all 
of those earlier experiments were carried out in a clas- 
sical regime at room temperature. Experiments in the 
quantum regime |7H_9| have long been anticipated, and 
the use of quantum-coherent mesoscopic conductors may 
lead to this goal. However, at their typical operation 
temperatures T below 1 K, the requirement to resolve 
tiny fluctuations on the energy scale fc^T becomes an in- 
creasingly challenging task. An interesting recent result 
in this direction has been the verification of exact rela- 
tions between current and current noise as functions of 
the source-drain voltage across an Aharonov-Bohm in- 
terferometer [10|. These relations are naturally derived 



from the FT [Eq. ([T])] as well as Onsager-Casimir and 
fluctuation-dissipation relations ll|, [iJl . 

As a step toward the direct test of Eq. ([T]) in the quan- 
tum regime, we verify the fluctuation theorem in single- 
electron tunneling 13l4l5[ at low temperatures, although 
our experiment is carried out in the regime of classical 
charge counting. We employ real-time detection of single- 
electron charging lallll in quantum dots (QDs). Mon- 
itoring the charge state of two QDs coupled both in se- 
ries and to source and drain electrodes (double quantum 
dot [DQD]) allows us to measure the direction- resolved 
charge flow through this device [ISJ and consequently the 
current probability distribution. A recent experiment 
following this line [19[ revealed the importance of the 
backaction '20| of the charge sensing device, which was 
a quantum point contact [21| (QPC). Backaction due to 
nonequilibrium QPC noise destroys microreversibility in 
the DQD and leads to an apparent temperature that sig- 
nificantly exceeds the bath temperature LlE,^,^]. This 
sensitivity to the measurement arises because the entropy 
fluctuations associated with single-electron tunneling are 
3 orders of magnitude smaller than those at room tem- 
perature. In order to avoid the spurious backaction, we 
employ an optimized sample design combining electron- 
beam and scanning-probe lithography |24| . It provides 



the high tunability and electronic stability required for 
the experiment while maintaining a good QPC-DQD 
coupling. We observe quantitative agreement between 
our data and theory in the near-equilibrium regime after 
including a small correction due to finite detector band- 
width [25i] , which is a technical rather than physical com- 
plication and requires no fitting parameter. In the regime 
far from equilibrium, our results depend on the system 
details; we find that the FT describes our data correctly 
in configurations where the DQD dynamics are those of 
a three-state Markovian system. 

Our measurements were performed in a "^He/^He di- 
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FIG. 1: (a) Atomic-force micrograph of the sample. Electrons 
can travel between the source and drain via the two quantum 
dots marked by disks. The conductance Gqpc of the quan- 
tum point contact serves to read out the charge state of the 
quantum dots, (b) By using the two gates Gl and G2, the 
charge state of the sample is controUably switched between 
empty (0), left (L), right (R) and doubly (2) occupied. These 
four states are visible as panels of distinct Gqpc in the plot, 
(c) Time dependence of Gqpc close to the charge degeneracy 
point marked by a dot in (b), displaying fluctuations between 
three levels: L, R, and 0. (d,e) Gqpc time segments show- 
ing examples of a direct transition from to L (d), and of a 
transition from to L via R (e). 



lution refrigerator on the sample shown in Fig. [Ija). 
The dark parts in the atomic-force micrograph corre- 
spond to the conductive (non-depleted) parts of a two- 
dimensional electron gas 34 nm below the surface of a 
GaAs/Alo.sGao.yAs heterostructure (sheet density ns = 
4.9 X 10^^ m~^, mobility ^ = 33m^V~-^s~^ as determined 
at r = 4.2 K). Confinement is achieved in one part with 
Ti/Au gates (the upper half of the image) biased with 
negative voltages. The thin vertical finger gates are only 
slightly biased in order to maintain a small tunneling cou- 
pling between the two QDs (white disks) and the source 
and drain leads. The horizontal lines are created by lo- 
cal anodic oxidation; they electrically separate the DQD 
from the charge detector QPC in the lower half of the im- 
age. The detector's conductance Gqpc is sensitive to the 
charge on the DQD and abruptly decreases if an electron 
is loaded to either of the two QDs. Figure [TJb) shows 
a color plot of the time-averaged QPC conductance as 
a function of the two gate voltages Vgi and Vg2, which 

26|. Four re- 



control the electron number on the DQD 
gions of stable charge are visible as regions of constant 
QPC conductance, up to a background that is linear in 
Vgi and Vg2- We estimate that each QD holds about 
80 electrons, but for ease of notation we consider only 



the excess electron number and denote the four relevant 
charge states of the DQD as empty (0), singly occupied 
(L, R), and doubly occupied (2). At the borders of these 
regions, thermal fluctuations of the charge occur. 

In particular, at the charge degeneracy point marked 
by a dot in the color plot, the QPC conductance rapidly 
fluctuates between three levels corresponding to the 
states 0, L, and R, as shown in Fig. [H^c). By counting 
the number of transitions between L and R, we can deter- 
mine the total number n of electrons that pass through 
the center barrier of the DQD during an acquisition time 
T. If T is large compared to the typical dwell time of an 
electron inside the DQD, the electron passing through the 
center barrier will typically reach one of the leads, where 
it equilibrates with the thermal bath at temperature T. 
The entropy change AS will then be equal for all charge 
fluctuations with equal n [9|, |l2|, |l5| . The change can be 
either positive or negative, depending on the direction of 
the charge flow with respect to the bias direction. As 
the dissipated heat neT^DQD is determined by the DQD 
source-drain voltage, AS* is given by neVoQu/T, and the 
FT for our system is 
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The QPC is biased with a voltage of 300 fiY and its 
conductance recorded using a room-temperature current- 
to-voltage converter and digitizer. The signal is then 
flltered by software at a bandwidth of 0.4 kHz, further 
resampled at 1.5 kHz, and stored for analysis. Two con- 
ductance thresholds were defined in the middle between 
neighboring conductance levels. In the algorithm to de- 
termine the chronological sequence of DQD states, we 
build in the requirement that at least three successive 
data points must lie in the same conductance interval for 
the DQD state to be accepted. The reason for this choice 
becomes clear when taking a closer look at the Gqpc time 
traces provided in Fig. [TJd) and (e) . The Gqpc signal in 
panel (d) exhibits a direct transition from level to level 
L. Because of to the finite rise time, two of the sampled 
data points happen to lie between thresholds 1 and 2. If 
these were to be assigned to level R, a false transition 
from R to L would be counted, reducing the net flow n 
in the corresponding time segment by 1. In comparison, 
the signal in panel (e) shows a short, yet clear, dwell time 
in the level R. 

Figure [2lja) shows an example of a Prin) distribution 
measured at an electron-bath temperature of 330 mK and 
with Vdqd = ^V. It is based on the counting analysis 
of 3000 Gqpc time segments, each with length r = 2s. 
The choice of t is such to minimize the combined error 
that originates both from the imperfect long-time limit 
[Sl, ll2, UM, uM (favoring large r) and from statistics (fa- 
voring small t). The distribution is symmetric about 
n = 0; i.e., there is no net charge flow, as expected in 
equilibrium. The amount of charge flow at zero DQD 
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FIG. 2: (a,b) Experimental probability distributions for the 
net electron number n transferred through the DQD during 
time r = 2s for two different DQD source-drain voltages 
Vdqd- (c,d,e) Comparison of experimental data with the- 
ory for three different bath temperatures. The data points 
correspond to the left-hand side of Eq. ((2]) and describe the 
probability ratio of forward (-|-n, entropy-producing) and 
backward (— n, entropy-consuming) processes for a given 
n. The solid lines mark the expected exponential behav- 
ior exp(nel/DQD/fcsT) for the two source-drain voltages O^V 
(dark blue, horizontal) and 20 /iV (red, inclined). If the fi- 
nite bandwidth of the detector is taken into account [23, [23] 
(dashed lines), experiment and theory agree within the sta- 
tistical uncertainty of the data. (Error bars indicate the esti- 
mated standard deviation. The gray bands around the dashed 
lines indicate the uncertainty in the finite-bandwicth correc- 
tion.) 



bias (compared to charge flow at finite DQD bias) is a 
sensitive measure for the strength of residual QPC back- 
action. Because of left-right asymmetries in the DQD 
such as a nonzero level detuning, QPC backaction gener- 
ically leads to nonequilibrium charge flow in a ratchet- 
type effect [iili^. 

When T^DQD is increased to 20 /iV, the distribution 
shifts towards positive n, as shown in panel (b), and, 
on average, an electron number of (n) = 2.5 is trans- 
ferred. Still, for some of the time segments, the charge 
flow is against the applied bias (n < 0), which results 
in a temporary decrease of the system entropy. Similar 
measurements were carried out at temperatures 500 niK 
and 700 mK fSlj. These are shown in Fig. [2ljc-e), where 
the data points are the logarithm of the left-hand side of 
Eq. ([2|), measured at Vdqd — fiV and Vdqd = 20 /zV, 
respectively. The expression ln[PT-(n)/PT-(— n)] follows 
the expected linear behavior close to the theoretical curve 
neVoQu/kBT (solid lines). 

In the nonzero-bias case, there is a systematic devi- 
ation of 20% to 30% in the slope. This can be under- 



stood by taking into account the limited bandwidth of 
the charge detection. A charge-switching event in the 
DQD is detected in the QPC only after a reaction time 
of l/Tdot, which in our case is determined both by the 
rise time of the measurement electronics and by the re- 
jection of short events built into the analysis algorithm. 
If the charge state switches back too fast, the event is 
missed. Following the ideas presented in Ref. [27|, Ut- 
sumi et al. [25| calculated the effect of the finite band- 
width. They found that, up to order l/Fdct, the finite 
bandwidth has the same effect as a prefactor aew < 1 to 
the term eVogo/kBT, just as we have observed in our 
experiment. The factor aew = ksT hi w* / cVdqd is ex- 
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where w = ejq){eVuQi)/kBT). Qualitatively, the effect 
can be understood as follows. At a nonzero bias, tran- 
sitions directed toward the drain occur with faster rates 
(rLo, Trl, Por) than those directed toward the source 
(Pol, Plr, Pro)- It is therefore more likely that the 
detector misses a charge flowing toward the drain than 
toward the source. In short, the ratio Pr{n) / Pr{—n) be- 
comes underestimated for n > 0. It is important to stress 
that, despite the similar phenomenology, the correction 
we employ here has a very different quality than the ef- 
fective temperature due to QPC backaction introduced in 
Ref. [l9j. QPC backaction implies excitation of degrees 
of freedom in the microscopic system of interest, whereas 
the effect of flnite bandwidth is a matter of the imperfect 
room-temperature electronics and leaves the microscopic 
system unperturbed. 

In the flnite-bandwidth model, Pdot plays the role of 
a mean inverse reaction time of the detector. We use a 
simulation to determine this parameter for our particu- 
lar detection scheme. A square pulse of duration Tp is 
flltered and resampled the same way as the experimental 
signal. Our value of Pdet — (0.59 ± 0.12) kHz is then de- 
flned as the inverse of the minimum Tp for which at least 
three sampled points reach 50 % of the square-pulse am- 
plitude. The dashed lines in Figs. [IJc-e) show the theo- 
retical expectation calculated with this value for Pdet and 
the experimentally determined Tij [l8| , and indeed agree 
much better with the experiment. The gray shaded areas 
indicate the uncertainty in the slope which is mainly de- 
termined by the uncertainty in Pdot- We emphasize that 
this analysis does not involve any free parameters. 

The DQD voltage of 20 ^V used in the temperature- 
dependence measurements is comparable with the ther- 
mal voltage ksT/e — 28...58//V, so the system is not 
too far from thermal equilibrium. The FT also applies 
far away from equilibrium, however. To test its pre- 
dictions in this regime, we have also performed bias- 
dependence measurements with Vdqd up to 120 /iV, i.e.. 
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FIG. 3: (a) The red data points show the Vdqd depencence 
of the left-hand side of Eq. Q, which is the ratio of entropy- 
consuming vs. entropy-producing cycles, with error bars that 
indicate its estimated standard deviation. The blue circles 
show the right-hand side, which is the average of the Boltz- 
mann factor among the entropy-consuming cycles. The FT 
[Eq. Q] is satisfied if the finite detector bandwidth is taken 
into account in the form of a correction factor aew to the 
exponent in Eq. Q (shown as crosses). The uncertainty in 
the finite-bandwidth correction is comparable with the error 
in Prin < 0)/Pt(h > 0) for all bias voltages, (b) Dwell-time 
distributions of the Gqpc signal in the states L, R, and 0, as 
extracted from the data set in (a) at the point Vdqd ~ 0/iV. 
(c,d) Same as (a) and (b), but measured at a slightly differ- 
ent gate-voltage configuration. In this case, the dwell times at 
zero bias (d) are nonexponentially distributed for the states L 
and R. The finite-bandwidth model [the crosses in (c)] is not 
valid for this case but plotted for comparison, (e) Diagram of 
the decay of a charge state L (black rectangle) that consists 
of two internal QD states LI and L2. Such a process can lead 
to nonexponentially distributed charge dwell times as shown 
in (d). 



about 4.3 X ksT/e at T = 330 mK. The data are shown 
in Fig. [3] and are based on the analysis of 2000 Gqpc 
time segments of length r = 2 s for each DQD voltage. 
For clarity and to reduce the statistical error for large 
voltages, we plot an integrated version of the FT relat- 
ing the total fractions of entropy-producing and entropy- 
consuming cycles, 
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In our measurements of the Vdqd dependence, we keep 
the gate voltages Vqi and Vg2 fixed. The choice of 
Vqi and Vg2 determines the level arrangement of the 
DQD with respect to the electrochemical potentials of 
the leads and is a priori not relevant for the validity 



of Eq. ^. In our measurements we observe a rather 
strong effect of this choice which is not fully understood. 
Thus, we plot in Figs. [SJa) and (c) the data of two 
representative measurements. The differing gate-voltage 
configurations in the two measurements result in differ- 
ing relative arrangements of the electrochemical poten- 
tials of source, drain, left, and right QD (denoted by 
Ms = eyDQD/2, /i_D = -eT^DQD/2, hl, and m_r, respec- 
tively). In measurement (a), we have fiL — P^r ~ 80/ieV 
and {^iL +M-r)/2 ~ 10/xeV, whereas in measurement (b), 
we have /xl — ^r ~ 55 /xeV and (//l -|- \J.n)l^ ~ —5 /.teV, 
as determined by finite-bias spectroscopy. 

The red data points in Fig. ^a) plot the left-hand 
side of Eq. (JH) . The quantity rapidly decreases with the 
voltage, as charge transfer against the bias occurs less 
and less frequently. Measurements at even higher DQD 
voltages are eventually limited by the necessary, expo- 
nentially increasing measurement time. The blue circles 
plot the right-hand side, calculated without the finite- 
bandwidth correction. Similar to the low-bias case, there 
is a systematic deviation. For the black crosses, the expo- 
nent in Eq. (|4]) is replaced with the bandwidth-corrected 
version aBwneVuQD / ksT , and we see that the observed 
deviation can entirely be attributed to this measurement 
issue. 

In the second measurement shown in Fig. [3jc), there 
is a larger discrepancy between the two sides of Eq. dU 
which goes beyond the effect of the finite detector band- 
width. We observe that this discrepancy coincides with 
the presence of nonexponential distributions of the ran- 
dom dwell times in the three DQD charge states. Figures 
^h) and (d) show the histograms of the signal dwell times 
in L, R, and measured at Vdqd = 0/iV for the two 
configurations. In configuration (a), the dwell times fol- 
low exponential distributions with a single lifetime. This 
is not the case in configuration (c), where the states L 
and R are not characterized by a single lifetime. This 
is indicative for the population of additional (excited) 
states in the DQD [30| • Considering the example of the 
left QD, a nonzero population pL2 of a long-lived excited 
state L2 means that the charge state L can no longer 
be identified with the QD ground state LI. The QPC, 
which monitors the dynamics of the charge states, sees 
a dwell time t in the state L distributed according to 
PLit^i exp(-i/iLi) +PL2tZl exp(-t/tL2), where ^li and 
ti,2 are the decay times of the two QD states into R and 
[cf. Fig. [SJe)]. In configuration (c), the QD energies 
are lower than in configuration (a), which makes pop- 
ulation of excited states more probable. Although the 
precise level arrangement cannot be reconstructed from 
our data, the measurement in Fig. ^c) demonstrates the 
sensitivity of such a test of the FT to the details of the 
DQD level structure. 

In conclusion, we have presented an extensive quanti- 
tative test of the fluctuation theorem for electron trans- 
port through a DQD, covering different temperatures and 



strong nonequilibrium transport conditions. Our results 
validate the theory in the near-equilibrium regime with 
a good accuracy. A remaining discrepancy is very well 
explained with a master-equation model of the finite- 
bandwidth detection. This agreement proves the use- 
fulness of this corrective approach in compensating for a 
slow detector. In the regime far from equilibrium, our 
results display a strong dependence on the internal DQD 
level structure controlled by gate voltages. In configu- 
rations where our system is well described as a three- 
state Markovian system, we observe a good agreement 
with theory, demonstrating the potential of the DQD as 
a model system for the study of nonequilibrium thermo- 
dynamics. Our results anticipate the test of the FT in 
quantum-coherent electron transport, which requires the 
measurement of thermal fluctuations on a sub-Kelvin en- 
ergy scale. 
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